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ABSTRACT 

Aims. We aim at studying the abundance gradients along the Galactic disk and their dependence upon several parameters: a threshold in the 
surface gas density regulating star formation, the star formation efficiency, the timescale for the formation of the thin disk and the total surface 
mass density of the stellar halo. We use chemical evolution models already tested on a large number of observational constraints. 
Methods. We test a model which considers a cosmological infall law. This law does not predict an inside-out disk formation, but it allows to 
well fit the properties of the solar vicinity. To check whether this model can reproduce the properties of the galactic disk, we study several cases. 
We drop the threshold in the surface gas density and assume a star formation efficiency varying with radius. We also test the same parameters 
in the two-infall model for the Galaxy. Finally, we perform some additional analysis on our simulations to test if the cosmological infall law 
can account for the inside-out formation of the disk. 

Results. We find that to reproduce at the same time the abundance, star formation rate and surface gas density gradients along the Galactic disk 
it is necessary to assume an inside-out formation for the disk. The threshold in the gas density is not necessary and the same effect could be 
reached by assuming a variable star formation efficiency. The derived new cosmological infall law contains a mild inside-out formation and is 
still not enough to reproduce the disk properties at best. 

Conclusions. A cosmologically derived infall law with an inside-out process for the disk formation and a variable star formation efficiency 
can indeed well reproduce all the properties of the disk. However, the cosmological model presented here does not have sufficient resolution to 
capture the requested inside-out formation for the disk. High resolution cosmological simulations should be performed to better capture such a 
behavior. 
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1. Introduction 

Understanding the formation and the evolution of the Milky 
Way is fundamental to improve the knowledge of the formation 
of spiral galaxies in general. In many models of the chemical 
evolution of the Milky Way gas infall has been invoked to ex- 
plain the formation of a galactic disk (e.g. Chiosi 1980; Lacey 
& Fall 1985; Matteucci & Francois 1989; Chiappini et al. 1997 
(hereafter C97); Boissier & Prantzos 2000, Portinari & Chiosi 
2000 among others). Originally, the gas infall was introduced 
as a possible solution to the G-dwarf problem (Pagel 1989). 
The infall law is clearly very important to determine the main 
characteristics of a galaxy. 

Our Galaxy has many observational constraints. Some of 
the most important ones are represented by the evolution of the 
abundances of the chemical elements. Recently, many chemical 
evolution models have been developed to explain the chemical 
composition of the solar vicinity (e.g. Henry et al. 2000; Liang 
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et al. 2001 ; Chiappini et al. 2003a, 2003b; Akerman et al. 2004; 
Frangois et al. 2004; Cescutti et al. 2006, 2007). Other impor- 
tant constraints are the abundance gradients of the elements 
along the disk of the Milky Way. Abundance gradients are ob- 
served in many galaxies with their metallicities decreasing out- 
ward from the galactic centers. They are strongly influenced 
by the star formation and by the accretion history as functions 
of galactocentric distance (see Matteucci & Francois 1989; 
Boissier & Prantzos 1999; Chiappini et al. 2001; Colavitti et 
al. 2008). In particular, Colavitti et al. (2008) (hereafter COS) 
adopted an infall law cosmologically motivated, in the sense 
that assumed an accretion law for baryons which is similar 
to the accretion law of dark matter, as derived by cosmologi- 
cal simulations, obtained with the public tree-code GADGET2 
(Springel 2005). This cosmological law is able to well repro- 
duce all the constraints in the solar vicinity, such as the star 
formation rate, the supernova (SN) la and SNII rate, the total 
amount of gas and stars and the behavior of several chemical 
elements (such as the a-elements, C, N and Fe). 
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In this paper we base our work on two chemical evolution 
models: the one of C97 and the one of COS. The COS model 
does not predict an inside-out formation of the thin disk as in 
the C97 model and while the COS model predicts good star 
formation rate and gas density distributions as functions of the 
galactocentric distance, it does predict a too flat slope for the 
oxygen gradient in the inner disk regions. 

The aim of this paper is to test the predictions of the COS 
model along the disk of the Milky Way by varying a number of 
important parameters in order to understand if there are alter- 
native solutions to the inside-out formation and the threshold 
in the gas density for star formation to reproduce abundance 
gradients. In particular, we test several cases in which we vary: 
i) the threshold gas density during the formation of the halo, 
thick and thin disk components, ii) the efficiency of star for- 
mation as a function of galactocentric distance, iii) the total 
surface mass density of the stellar halo and iv) the timescale 
for the formation of the thin disk. We do not test the eff'ects of a 
variable initial mass function (IMF) along the disk, since it has 
already been done, and it has been concluded that it cannot re- 
produce at the same time the abundance gradients, the gas and 
star formation rate distributions (e.g. Carigi 1996; Chiappini et 
al. 2000). Radial flows can also play a role in the formation 
of abundance gradients (see Portinari & Chiosi 200, and refer- 
ences therein). In particular, radial flows can produce gradients 
but mainly in the outer regions of the disk while the inner gra- 
dients are too flat. This effect is similar to that of the threshold 
in the gas density without an inside-out formation and therefore 
we do not take it into account here. What is more important is 
the eff'ect produced by a bar in the bulge which can explain the 
peak of the gas observed at 4-5 Kpc from the Galactic center 
(Portinari & Chiosi 2000), contrarily to models which do not 
include the bar. However, we are not interested in this problem 
here but simply to the gradients and gas and star formation rate 
distributions for a galactocentric distance Re > 4 kpc. 

We compare our model predictions to observational data 
relative to abundance gradients along the disk: those of 
Andrievsky et al. (2002a-c, 2004) and Luck et al. (2003). They 
measured the abundances of all the selected elements (O, N and 
Fe) in a sample of 130 Galactic Cepheids found in the galacto- 
centric distance range from 5 to 17 kpc. Moreover, we compare 
our results with the observational data collected by Simpson 
et al. (1995), Afilerbach et al. (1997) and Gummersbach et al. 
(199S) for the N gradient from Hll regions and with those col- 
lected by Chen, Hou & Wang (2003), Twarog et al. (1997) and 
Carraro et al. (1998) for the Fe gradient from open clusters. We 
also compare our model results with the star formation rate and 
the surface gas density as functions of the galactocentric dis- 
tance by adopting the data collected by Rana (1991). He took 
the average HI distribution from Blitz et al. (19S3), Gordon 
& Burton (1976) and Garwood & Dickey 19S9 and two dis- 
tinct H2 distributions: one from Bronfman (1986) and the other 
from Robinson et al. (19SS), Wouterloot et al. (1990) and Digel 
et al. (1990). Moreover he took the radial dependences of the 
present surface density of the star formation rate from Rana & 
Wilkinson (19S6). 

The paper is organized as follows: in section 2 we show the 
observational data and nucleosynthesis prescriptions adopted. 



Section 3 presents a brief description of the model by C97. 
In section 4 we describe the COS model and the cosmological 
simulations. Section 5 describes the results obtained, compar- 
ing the models predictions with the observed properties. FinaUy 
section 6 presents the conclusions. 

2. Observational data and nucleosynthesis 
prescriptions 

In this work we use the data by Luck et al. (2003) and 
Andrievsky et al. (2002a-c, 2004) for all the studied elemen- 
tal gradients. These data have been derived for a large sam- 
ple of Galactic Cepheids. These variable stars have a distinct 
role in the determination of radial abundance gradients for sev- 
eral reasons. First, they are usually bright enough that they can 
be observed at large distances, providing accurate abundances; 
second, their distances are generally well determined, as these 
objects are often used as distance cahbrators; third, their ages 
are also well determined. They generally have ages close to a 
few hundred million years. So we can assume that they are rep- 
resentative of the present day gradients. 

Here we also use the data by Gummersbach et al. (1998) (B 
stars), Simpson et al. (1995) and Afilerbach et al. (1997) from 
HII regions, and by Chen, Hou & Wang (2003), Carraro et al. 
(1998) and Twarog et al. (1997) from open clusters. 

One of the most important ingredients for chemical evolu- 
tion models is represented by the nucleosynthesis prescriptions 
and, consequently, by the stellar yields. 

The single stars in the mass range O.S < M < S M© 
(low and intermediate-mass stars) contribute to the Galactic en- 
richment through planetary nebula ejection and quiescent mass 
loss. They enrich the interstellar medium mainly in He, C, N 
and heavy s-process elements (e.g. Cescutti et al. 2006). We 
adopt here the stellar yields for low and intermediate mass 
stars of van den Hoek & Groenewegen (1997) computed as 
functions of stellar metallicity, their case with variable mass 
loss. Concerning N production, the yields of van den Hoek & 
Groenewegen (1997) from AGB stars predict both primary and 
secondary N. For N produced in massive stars we assume it 
secondary, although there is some indication that massive stars 
might produce primary nitrogen (e.g. Meynet & Maeder 2002) 
A small fraction of these stars are also the progenitors of Type 
la SNe if they are in binary systems: they originate from carbon 
deflagration of C-O white dwarfs. We adopt in this paper the 
single-degenerate progenitor scenario (Whelan & Iben, 1973; 
Han & Podsiadlowski 2004). Type la SNe contribute a substan- 
tial amount of Fe (~ 0.6 Mq per event) and Fe-peak elements 
as well as non negligible quantities of Si and S. They also pro- 
duce other elements, such as O, C, Ne, Ca, Mg and Ni, but in 
very small amounts compared to Type 11 SNe. We assume the 
stellar yields for Type la SNe from Iwamoto et al. (1999). 

Massive stars (8 M© < M < 100 M©) are the progenitor 
of core-collapse SNe which can be either Type 11 SNe or Type 
Ib/c SNe. These latter can arise either from binary systems or 
Wolf-Rayet stars, whereas Type II SNe originate from the mas- 
sive stars in the lower mass range (Calura & Matteucci 2006). 
Type 11 SNe mainly produce the so called a-elements, such as 
O, Mg, Ne, Ca, S, Si and Ti, but also some Fe and Fe-peak 
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elements although in smaller amounts than Type la SNe. We 
adopt here the stellar yields for massive stars by Woosley & 
Weaver (1995) with the suggested modifications of Francois et 
al. (2004). However, the most important modifications concern 
some Fe-peak elements, except Fe itself, whereas for the a- 
elements, with the exception of Mg which has been increased 
relative to the original yields, the yields are substantially un- 
modified relative to Woosley & Weaver's ones. 

The reference solar abundances are those by Asplund et al. 
(2005). 

3. The model by Chiappini et al. (1997) 

The model by C97 was the first in which two main infall 
episodes for the formation of the Galactic components were 
suggested. In particular, they assumed that the first infall 
episode was responsible for the formation of the halo and thick- 
disk stars that originated from a fast dissipative collapse. The 
second infall episode formed the thin-disk component, with a 
timescale much longer than that of the thick-disk formation. 
The timescale for the formation of the halo-thick disk is 0.8 
Gyr, while the timescale for the thin disk is 7 Gyr in the solar 
vicinity. The authors included in the model also a threshold in 
the gas density, below which the star formation process stops. 
The existence of such a threshold value is suggested by obser- 
vations relative to the star formation in external disk galaxies 
(Kennicutt 1998, but see Boissier et al. 2006). The physical 
reason for a threshold in the star formation is related to the 
gravitational stability, according to which, below a critical den- 
sity, the gas is stable against density condensations and, con- 
sequently, the star formation is suppressed. In the two-infaU 
model the halo-thick disk and the thin disk evolutions occur at 
different rates and they are independent, mostly as a result of 
difl^erent accretion rates. With these precise prescriptions it is 
possible to reproduce the majority of the observed properties 
of the Milky Way and this shows how important is the choice 
of the accretion law for the gas coupled with the star formation 
rate in the Galaxy evolution. Some of the most important ob- 
servational constraints are represented by the various relations 
between the abundances of metals (C, N, a-elements, iron peak 
elements) as functions of the [Fe/H] abundance and by the G- 
dwarf metalhcity distribution. In C97 the Galactic disk is ap- 
proximated by a series of concentric annuli, 2 kpc wide, with- 
out exchange of matter between them. In this model the thin 
disk forms "inside-out", in the sense that it assumes that the 
timescale for the disk formation increases with galactocentric 
distance. This choice was dictated by the necessity of repro- 
ducing the abundance gradients along the Galactic disk. The 
SFR is a Schmidt (1955) law with a dependence on the surface 
gas density {k = 1.5, see Keimicutt 1998) and also on the total 
surface mass density (see Dopita & Ryder 1994). In particular, 
the SFR is based on the law originally suggested by Talbot & 
Arnett (1975) and then adopted by Chiosi (1980): 

/2(r,0V'(r,/)V*""^ , ,k 

where t is the time and the constant v is a sort of efliciency 
of the star formation process and is expressed in Gyr~^: in 



particular, v = 2 Gyr~^ for the halo and 1 Gyr~^ for the disk 
(t > I Gyr). The total surface mass density is represented by 
£(r, t), whereas 2(rQ, t) is the total surface mass density at the 
solar position, assumed to be = 8 kpc (Reid 1993). The 
quantity 'Lgasi'', f) represents the surface gas density and t rep- 
resents the time. These choices of values for the parameters 
allow the model to fit very well the observational constraints, 
in particular in the solar vicinity. A threshold gas density for 
the star formation in the thin disk of JMqPc^^ is adopted in 
all the models presented here. The IMF is that of Scalo (1986) 
normalized over a mass range of 0.1-100 M© and it is assumed 
to be constant in space and time. In this paper we will adopt 
also a simple Schmidt law where the SFR depends only on the 
surface gas density: 

ifrir,t)=v-Lgasir,tf (2) 
where k = 1 .5 and v is proportional to 1/R. 

4. The model of Colavitti et al. (2008) 

The model of C08 is basically the one of C97 but with an ac- 
cretion law derived from cosmological simulations. In particu- 
lar, we run a dark matter-only cosmological simulation, using 
the public tree-code GADGET2 (Springel 2005), in order to 
produce and study dark matter halos in which spiral galaxies 
can form. Our simulated box has a side of 24 ' Mpc. We 
used 256^ particles. We adopted the standard cosmological pa- 
rameters from WMAP 3-years (Spergel et al. 2007), namely 
Qo = 0.275, = 0.725 and il^ = 0.041. Each DM parti- 
cle has a mass equal to 6.289 • lO^/i'^M© and the Plummer- 
equivalent softening length is set to 3.75 comoving kpc to 
redshift z - 2 and to 1.25 /z ' physical kpc since z = 2. We 
use the pubUc package GRAFTC (Bertschinger 1995) to set up 
our initial conditions. The simulation started at redshift z = 20 
and 28 outputs were produced. We have chosen to use a quite 
large spread in the redshifts at the beginning, while in the last 
part of the simulation, where a small change in the redshift cor- 
responds to a large change in time, the redshifts are closer. We 
checked that the final mass function of DM halos and the power 
spectrum are in agreement with theoretical expectations. 

We identified DM halos at redshift z = using a stan- 
dard friend-of-friends algorithm, with a linking length / = 0.17 
mean (comoving) interparticle distance. After that, we deter- 
mined the virial mass and radius for each DM halo, using the 
center of mass of the F-o-F group as the halo center. Here we 
define the virial radius as the radius of the sphere within which 
the matter density contrast is6 ^ 100 times the critical density, 
with S given by the cosmological parameter as in Navarro & 
Steinmetz (2000). 

We then built the mass accretion history of our halos. To 
achieve this goal, we analysed the 28 outputs obtained from 
redshift z = 9.0 to z = 0. We identified all DM halos in each 
snapshot using the procedure described above, except for the 
fact that we used the redshift-dependent density contrast given 
by Bryan & Norman (1997) to define the virial radius as a func- 
tion of z. At any output Zi+i , we found all the progenitors of our 
halos at redshift z,. We defined a halo at redshift Zi+i to be a 
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progenitor of one at z, if at least 50% of its particles belong to 
the candidate offspring (see e.g. Kauffmann 2001, Springel et 
al. 2001 for a discussion of this threshold). The mass accretion 
history is defined as the mass of the main progenitor of the halo 
as a function of redshift. Having the mass accretion histories, 
we were able to identify the redshift of formation (defined as 
the epoch at which half of the mass of the forming halos were 
accreted) and the redshift at which each halo experienced its 
last major merger (defined as an increase of at least 25% of its 
mass with respect to the mass of its main progenitor at the pre- 
vious redshift). To identify the DM halos that can host a spiral 
galaxy similar to the MW we used selection criteria based on 
four different characteristics of the halos: 

- mass between 5 ■ IO^Mq and5 ■ IO'^Mq; 

- spin parameter A > 0.04; 

- redshift of last major merger (Zmm) larger than z = 2.5; 

- redshift of formation (zf) larger than z - 1 .0. 

In particular, in this paper we use our best cosmological 
halo which has these characteristics: 

- M = 90.26 • IO'^Mq; 

- 0.045; 

~~ Zmm~ 5.00, 

- 1.50 <z/< 1.75. 

Assuming that the baryonic matter follows the same accre- 
tion pattern as the dark matter, and that it represents 19% (the 
cosmological baryon fraction) of all the infalling matter, we 
obtained a final baryonic mass for the Galaxy of 1.7 ■ lO^M©. 
We assumed that the derived infall law has the same functional 
form for the whole Milky Way, but that the normalization con- 
stant is different for different Galactic regions. In other words, 
the normalization constants were obtained by reproducing the 
current total surface mass density at any specific galactocentric 
distance. 

5. Results 

In this work we tested several different models in order to un- 
derstand if it is possible to well reproduce all the observed val- 
ues just acting on four variables. In particular, in COS (Model 
COS) we saw that our best cosmological model was not able 
to reproduce the [O/H] ratio in the inner Galactic disk. So in 
this paper we try to find the better combination of variables to 
reproduce the oxygen gradient. We also test these variables in 
the model by C97. These are the considered variables: 

- Threshold. In the original model by C97 the threshold in 
the surface gas density was equal to 4 Mq pc^^ during the 
formation of the halo and to 7 pc"^ during the for- 
mation of the thin disk. When the surface gas density was 
lower than these values the SFR stopped. In this paper we 
switch off the threshold in all the models. 

- SF efficiency (v). In C97 and in COS the star formation ef- 
ficiency was assumed to be constant along the disk. Here 
we use, for some models, a SF efficiency which is a func- 
tion of radius. This approach is already known (see Boissier 
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Fig. 1. Cosmological infall laws at 4 (blue dotted line), S (red 
solid line) and 14 kpc (green dashed Une). 



& Prantzos 1999) and such a dependence has been pro- 
posed on the basis of large-scale instabilities in rotating 
discs. Even in Carigi (1996) the author used a model in 
which the star formation efficiency was variable with radius 
(v(r) oc r"' , as akeady seen in Prantzos & Aubert 1995). We 
adopt several ad hoc variations of the parameter v as a func- 
tion of radius plus the law adopted in Boissier & Prantzos 
(1999). All the adopted laws are shown in Table 1. 
Ximescale for the formation of the disk (r). In C97 the 
Galactic disk forms inside-out. This means that the in- 
ner parts of the galactic disk formed first than the outer 
ones. This type of formation was introduced to explain the 
Galactic gradients. Even in Carigi (1996) the author used a 
model in which the infall timescale for the disk was vari- 
able with radius. The model by COS, instead, did not use 
an inside-out scenario since the cosmologically derived in- 
fall law does not contain such process. So in this case the 
parameter r cannot be changed. In Figure 1 we show the 
cosmological infall law of COS at different radii. Clearly 
the amount of gas falling at each radii is different, but the 
shape of the law is maintained. 

In some of the models we use a constant timescale while in 
others we use the prescription used in C97, namely: 



T= 1.03 r- 1.27 [Gyr]. 



(3) 



where r is the galactocentric radius in kpc. Other authors 
have also studied an inside-out formation of the disk by 
assuming a linear dependence of t on r (e.g. Matteucci & 
Frangois 19S9, Carigi 1996, Boissier & Prantzos 1999). 
Surface halo mass density (I.haio)- In C97 the projected 
halo mass density was a constant parameter. Here we use 
both a halo density which is a function of radius and a con- 



E. Colavitti et al.: Chemical evolution along radius 



5 



Slant one, as in Chiappini et al. (2001). In particular, it fol- 
lows this equation: 

I^haio ^ — [Me pc-^]. (4) 
r 

where r is the radius in kpc. The value of 136 is adopted 
in order to have a surface halo mass density equal to 17 
Mq pc^^ at the solar neighbourhood. The halo density is 
important in determining the slope of the gradient at large 
galactocentric distances, where the disk density is very low 
(see Chiappini et al. 2001). 

In Table 1 we present the models and all the parameters 
which characterize them. The second column indicates the 
presence or the absence of a threshold, the third column the 
type of variation of the SF efficiency (the values of v for Model 
3 are for 2.5, 7.5 and 12.5 kpc), the fourth shows if the galaxy 
forms inside-out or not and the fifth if the halo density is con- 
stant or if it is a function of radius. In the sixth column we show 
the type of infall law used and in the seventh column the SFR. 

Model COS is the best cosmological model of Colavitti et 
al. (2008). It had the threshold only during the formation of the 
disk and a constant star formation efficiency v. As explained 
above, with this model the authors were not able to reproduce 
the observed values for the [O/H] in the inner part of the disk. 

Model C97 is the original model of C97, i.e. the two-infall 
model. It has a threshold, both during the formation of the halo 
and thick-disk component and during the formation of the thin 
disk. It has a star formation efficiency constant along radius 
and equal to v = 1.0 Gyr an inside-out prescription and a 
constant surface halo density along radius. 

5.1. Results for models COS, C97, 1 and 2 

In Figures 12] and [3] the results for Model C97, compared with 
those predicted by Model COS and with observations, can be 
seen. Figure |2] upper panel, presents the SFR normalized at 
the solar neighbourhood value. The data are from Rana (1991) 
(green dotted line) while the red dashed line represents Model 
COS. The black solid line represents Model C97. It can be seen 
that both Model C97 and Model COS have no SF in the outer 
part of the disk, both having a threshold in the surface gas den- 
sity below which the SF stops. The bottom panel of Figure |2] 
shows the current surface gas density along the radius. The 
data are, once again, from Rana (1991). In this case Model C97 
and Model COS produce very similar results. Moreover both the 
models reproduce very well the lower values of Rana (1991). 

Figure |3] shows the evolution of [O/H], [N/H] and [Fe/H] 
along the radius for Models C97 and COS. The model results 
are normalized to the Asplund et al. (2005) solar abundances in 
all the models. The data are from a compilation by Cescutti et 
al. (2007) (blue dots, Cepheids stars), from Gummersbach et al. 
(199S) (red squares, B stars), from Simpson et al. (1995) and 
Afflerbach et al. (1997) (green triangles, HII regions) and from 
Chen, Hou & Wang (2003), Carraro et al. (199S) and Twarog 
et al. (1997) (cyan crosses, open clusters). The black squares 
are the mean values inside each bin only for the compilation by 
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Fig. 2. S FR/S FRq and surface gas density vs radius for Model 
C97 (black solid line) and for Model COS (red dashed Hne). The 
data are from Rana (1991) (green dotted lines, representing the 
lower and the maximum values for the gas). 



Cescutti et al. (2007) and the error bars are the standard devia- 
tions. It can be seen that Model C97 reproduces quite well the 
data, even in the inner part of the disk, since the slope is more 
pronunced than that of Model COS. As already said. Model COS 
is too flat in the inner part of the disk. For this reason we test 
several parameters in order to better reproduce the inner abun- 
dance gradients. In the outer parts of the galactic disk the slope 
of Model COS is more pronunced than that of Model C97. 

Figures |4] and |5] show the results for Model 1 , compared 
with those predicted by Model COS and with observations. 
Model 1 is equal to Model COS, except for the lack of the 
threshold. We can see that Model 1 reproduces very well the 
SFR/SFRq, much better flian Model COS. However, Model 1 
can not reproduce the surface gas density in the outer parts of 
the disk, since the obtained values are too low respect to the 
observed values. 

Figure |5] shows the abundance gradients for Model 1, com- 
pared with the results of Model COS. The data are the same 
of Figure [3] In this case the lack of a threshold in the surface 
gas density makes the abundance gradients completely flat. 
Therefore, simply removing the threshold from Model COS, it 
is not possible to reproduce the chemical abundances along the 
galactic disk, neither in the inner parts nor in the outers. So in 
the next model (Model 2) we also use a variable star formation 
efficiency along radius. 

In Figures |6] and |2] the results obtained for Model 2, com- 
pared with those predicted by Model COS and with observa- 
tions, can be seen. Model 2 has the same infall law of Model 
COS but in this case there is no threshold during the formation 
of the galactic disk. Moreover, in Model 2 we use a star forma- 
tion efficiency which is a function of radius. 
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Table 1. Models parameters. The second column indicates the presence or the absence of a threshold, the third column the type 
of SF efficiency, the fourth if the galaxy forms inside-out and the fifth if the halo density is constant or if it is a function of radius. 
In the sixth column we show the type of infall law used and in the seventh column the SFR. The values of v for Model 3 are for 
2.5, 7.5 and 12.5 kpc. 
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Fig. 3. [O/H], [N/H] and [Fe/H] vs radius for Model C97 (black 
solid line) and for Model COS (red dashed line). The data 
are from a compilation by Cescutti et al. (2007) (blue dots, 
Cepheids), from Gummersbach et al. (1998) (red squares, B 
stars), from Simpson et al. (1995) and Afflerbach et al. (1997) 
(green triangles, HII regions) and from Chen, Hou & Wang 
(2003), Carraro et al. (1998) and Twarog et al. (1997) (cyan 
crosses, open clusters). The black squares are the mean val- 
ues inside each bin only for the compilation by Cescutti et al. 
(2007) and the error bars are the standard deviations. 

In Figure |6] upper panel, it can be seen that Model 2 is in 
better agreement with the data then Model C08, especially in 
the outer part of the disk. Model C08 has aS FR/S FRq equal to 
zero in the outer parts of the disk, and it is due to the presence 
of the threshold. From the bottom panel of Figure |6] we can 
see that Model C08 is in quite good agreement with the data, 
whereas Model 2 predicts a surface gas density at variance with 
the observations: it predicts a too low surface gas density in the 
inner disk because of the high SF efficiency. 



Fig. 4. S FR/S FRq and surface gas density vs radius for Model 
1 (black solid line) and for Model C08 (red dashed line). The 
data are the same of Figure |2] 



Figure El shows the evolution of [O/H], [N/H] and [Fe/H] 
along the radius for Models C08 and 2. Using no threshold and 
a star formation efficiency variable with radius again we cannot 
reproduce, as we hoped, the gradients in the inner parts of the 
disk. We can only increase the slope of the gradients, relative 
to Model C08 in the outer parts. 

Since it is not possible to well reproduce at the same time 
the surface gas density and the abundance gradients along the 
galactic disk using the cosmologically derived infall law by 
C08 without the threshold or without the threshold and adopt- 
ing a SF efficiency variable along radius, we now try a new 
approch. As we said in section 4, model C08 assumed that the 
derived infall law has the same functional form for the whole 
Milky Way, but that the normalization constant is different for 
different Galactic regions. 
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Fig. 5. [O/H], [N/H] and [Fe/H] vs radius for Model 1 (black 
solid line) and for Model COS (red dashed line). The data are 
the same of Figure |3] 



Fig. 7. [O/H], [N/H] and [Fe/H] vs radius for Model 2 (black 
solid line) and for Model COS (red dashed line). The data are 
the same of Figure[3] 



3 - 
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Radius [kpc] 

Fig. 6. S FR/S FRq and surface gas density vs radius for Model 
2 (black solid line) and for Model COS (red dashed line). The 
data are the same of Figure|2] 

5.2. Results for the cosmological infall law with 
inside-out 

We then tried to derive infall laws having different functional 
form along the disk. To do this we divided the galactic disk in 
three zones. The first includes the central 5 kpc of the disk, the 
second is a radial shell starting at 5 kpc and ending at 10 kpc 
whereas the third shell has inner and outer radii of 10 and 15 
kpc. Having the coordinates of each DM particle of the main 
progenitor at every redshift, we calculated the total mass in- 



cluded in every shell, deriving three different cosmological in- 
fall laws for the three zones. It is well-known (Helmi et al. 2003 
and references therein) that the build-up of a Galaxy-sized DM 
halo proceeds in an inside-out fashion, with the mass of the in- 
ner part of the halo being in place at high redshift, while its 
outer part still accretes mass up to low redshift. We wanted to 
determine if a cosmological, standard resolution simulation can 
still capture such a behaviour. This in order to verify if it is pos- 
sible to use such a simulation to build a cosmological inside- 
out barionic accretion law, without resorting to high-resolution 
resimulation of single halos with DM only or with more com- 
ples physics. It is also known that re-simulation of single ha- 
los, with gas physics, a star formation prescription and some 
effective form of SN energy feedback can reproduce an inside- 
out formation of disk galaxies (Governato et al. 2007, Abadi 
et al. 2003a,b). Such numerical simulations are however com- 
putationally very expensive and the properties of the simulated 
galaxy somehow depend on the sub-grid prescription for the 
various astrophysical processes. Our aim here is to see whether 
we can obtain similar results using a much simpler numerical 
approach. 

Then we tested these laws in the chemical evolution model, 
using no threshold and a star formation efficiency which 
changes with radius. In Figures[8]and|9]the results can be seen. 

Figure [8] shows that Model 3 reproduces very well the ob- 
served SFR, especially between 9 and 12.5 kpc. The slope of 
the curve of observable data is very similar to that of Model 
3. Only at 2.5 kpc the predicted SFR is slightly larger then the 
observed. However the improvement respect to model COS is 
clear and it is mainly due to the lack of a threshold in the su- 
perficial gas density. In the bottom panel it can be seen that 
Model 3 has a too low "Egas respect to the observable data and 
model COS. However it is interesting to note that the behavior 
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Fig. 8. S FR/S FRq and surface gas density vs radius for Model 
3 (black points) and for Model COS (red dashed line). The data 
are the same of Figure |2l 



Fig. 9. [O/H], [N/H] and [Fe/H] vs radius for Model 3 (black 
solid line) and for Model COS (red dashed line). The data are 
the same of Figure[3] 



of the curve is quite similar to that of the observable data since 
the increasing of the superficial gas density observed by Rana 
(1991) in the outer part of the disk is also predicted by Model 
3. Moreover the value at 2.5 kpc is equal to that given by the 
observable data. 

Figure|9]shows the obtained [O/H], [N/H] and [Fe/H] along 
the radius for Model 3. It can be seen that in the case of O and 
N the results are quite similar to those of model COS, even if 
the slope produced by Model 3 is more pronunced. The slope is 
even more evident in the case of [Fe/H], especially for the outer 
part of the disk. However, not even Model 3 can well reproduce 
the slope in the inner part of the Galaxy, in spite of the fact that 
the accretion law contains an inside-out eff'ect. This is due to 
the mild inside-out effect predicted by the infall law which can 
be attributed to the relatively low resolution that our cosmolog- 
ical simulation can achieve on a single halo, especially at high 
redshifts. 

5.3. Changing the parameters in the C97 models 

After this study on the cosmological model by COS, we consid- 
ered the model by C97, trying to find the best combination of 
parameters to well reproduce the S FR/S FRq, the surface gas 
density and the abundance gradients along the Galactic disk, 
as done with Model COS. In this case we studied the behav- 
ior of the parameters in four different models, comparing every 
model with Model COS. 

Figures[TO]and[TT]show the results for Model 4. This Model 
has no threshold and the SF efficiency v, t and Sfta/o are con- 
stants. In Figure [To] it can be seen that the predicted star for- 
mation rate of Model 4 agrees very well with the data by Rana 
(1991). However l-gas is too low in the outer parts of the disk. 




4 6 8 10 12 14 



Radius [Kpc] 

Fig. 10. S FR/S FRq and surface gas density vs radius for 
Model 4 (black soUd line) and for Model COS (red dashed fine). 
The data are the same of Figure |2] 



even if in the inner parts it agrees quite well with the observed 
values. 

Figure [TT] shows the [O/H], [N/H] and [Fe/H] along radius 
for Model 4. We can see that in the inner disk the abundance 
gradients are flat, while in the outer parts they even increase, in 
contrast with the observed data. Therefore, it is not possible to 
reproduce the abundance gradients without the threshold and 
when V, t and l,i,aio are constants. This is due to the fact that 
with a constant t the disk evolution is too quick to allow the 
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Fig, 11. [O/H], [N/H] and [Fe/H] vs radius for Model 4 (black 
solid line) and for Model COS (red dashed line). The data are 
the same of Figure [3] It can be seen that in the outer part of 
the disk the gradients increase. In fact, being the halo surface 
mass density constant along the galactic disk and since there is 
no threshold in the surface gas density, the star formation in the 
halo overcomes that in the outer disk. 



formation of abundance gradients, that the star formation never 
stops (absence of the threshold) and that the halo surface mass 
density is constant with galactocentric distance. In this case, in 
fact, the star formation in the halo overcomes that in the outer 
disk and produces higher abundances than in the more realistic 
case of a fading surface halo mass density. 

FiguresfT2]and[T3]show Model 5, compared to Model COS. 
Model 5 has no threshold (neither in the halo nor in the disk), 
a constant star formation efficiency, a variable r along the disk 
and a constant surface halo density. From Figure [12] it can be 
seen that the star formation rate predicted by Model 5 is in quite 
good agreement with the data by Rana (1991), especially for 
the outer parts of the galaxy. However the amount of gas is too 
low, even if compared with the lower values by Rana (1991). 

Figure [13] shows the abundance gradients along the galac- 
tic disk for Model 5. In the inner part of the disk the [O/H], 
[N/H] and [Fe/H] ratios are not flat as in Model COS, reproduc- 
ing better the observed data. Nevertheless, beyond 10 kpc the 
gradients become positive at variance with the observations. So 
we can say that a model without threshold and with a variable 
T cannot well reproduce the data, especially in the outer disk. 

Figures [14] and [15] show the S FR/S FRq, "Egas and the 
[O/H], [N/H] and [Fe/H] gradients for Model 6. This model 
has no threshold, a star formation efficiency which changes 
strongly with radius, an inside-out prescription and a surface 
halo density which is a function of radius. From Figure [14] it 
can be seen that the SFR is in good agreement with the data, 
especially in the outer parts of the disk. On the other hand, the 
surface gas density increases along radius, at variance with the 



Fig. 12. S FR/S FRq and surface gas density vs radius for 
Model 5 (black sohd line) and for Model COS (red dashed hne). 
The data are the same of Figure [2] 
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Fig. 13. [O/H], [N/H] and [Fe/H] vs radius for Model 5 (black 
solid line) and for Model COS (red dashed line). The data are 
the same of Figure [3] This model shows that the presence of a 
variable t along the galactic disk steepens the inner gradients, 
without however reducing the positive slope in the outer disk, 
as seen in Fig.fTTI 



observations which show that there is more gas in the inner 
parts of the disk. 

In Figure [15] we can see that the abundance gradients have 
a slope which is quite similar to the slope of the data, mainly 
for the [Fe/H] ratio. Even in the inner parts the gradients are 
in quite good agreement with the data. However, we cannot 
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Fig. 14. S FR/S FRq and surface gas density vs radius for 
Model 6 (black solid line) and for Model COS (red dashed line). 
The data are the same of Figure |2] 
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Fig. 16. S FR/S FRq and surface gas density vs radius for 
Model 7 (black solid line) and for Model COS (red dashed hne). 
The data are the same of Figure |2] 
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Fig. 15. [O/H], [N/H] and [Fe/H] vs radius for Model 6 (black 
sohd Une) and for Model COS (red dashed line). The data are 
the same of Figured 



say that this model can reproduce all the observations, since it 
cannot reproduce the correct behavior of I,gas- 

The last model we present is Model 7. Its results are shown 
in Figures [T6l and [TtI This model has no threshold, a star for- 
mation efficiency which follows that one given by Boissier & 
Prantzos (1999) and no inside-out prescription. In this model 
we assume only one episode of accretion where halo forms first 
and then the disk follows. It differs from all the other models 
also because it adopts a simple Schmidt law for star formation. 



like in eq. (2). It can be seen that the SFR agrees well with the 
data by Rana (1991), while the surface gas density is too low 
in the inner part of the Galactic disk. This is again due to a too 
high SF efficiency in the inner disk. 

From Figure [17] we can see that Model 7 predicts abun- 
dance gradients which are in good agreement with the ob- 
served data, especially for the outer parts of the Galactic disk. 
Therefore we can say that this is the only case with no thresh- 
old and no inside-out prescription where it is possible to well 
reproduce abundance gradients, although the surface gas den- 
sity is too low in the inner part of the Galaxy. 

6. Conclusions 

Starting from the model by Colavitti et al. (200S), adopting a 
cosmologically motivated infall law resembling the one of C97 
for the solar neighbourhood but not predicting an inside-out 
formation of the Galactic disk, and in which the authors were 
not able to reproduce the [O/H] gradient in the inner parts of 
the disk, and from the model by C97, we tested several models 
in order to find the best values of four parameters to best re- 
produce the SFR/SFRq, Ygas and the [O/H], [N/H] and [Fe/H] 
gradients along the Galactic disk. These four parameters are the 
threshold gas density for the star formation, the star formation 
efficiency v, the inside-out prescription T(r) (i.e. the variation 
of the timescale for the formation of the Galactic disk) and the 
surface halo mass density. 

In table 2 we summarize model successes and failures 
for abundance trends, surface gas densities and star formation 
rates. 

Our conclusions can be summarized as follows: 

- We found that it is impossible to fit all the disk constraints 
at the same time without assuming an inside-out formation 
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Model 



Successes 



Failures 



good agreement for the S FR/S FRq 



"Lgas not reproduced in the outer Galactic disk 
abundance gradients completely flat 



good agreement for the S FR/S FRq 
increased slope of the gradients in the outer disk 



too low Zga, in the inner disk 
abundance gradients too flat in the inner disk 



good agreement for the S FR/S FRq 
Ejoi increases in the outer part of the disk 
slope of abundance gradients more pronunced in the outer disk 



too low Ejo, 
abundance gradients too flat in the inner disk 



good agreement for the S FR/S FR^ 



S^os not reproduced in the outer Galactic disk 
abundance gradients increase along the galactic disk 



good agreement for the S FR/S FRq 
slope of abundance gradients more pronunced in the inner disk 



Ego, not reproduced in the outer Galactic disk 
abundance gradients increase in the outer Galactic disk 



good agreement for the S FR/S FRq 
abundance gradients in good agreement with the observed data 



l,gas increases along radius 



good agreement for the S FR/S FRq 
abundance gradients in good agreement with the observed data 



Sgfl, too low in the inner part of the Galactic disk 



0.5 - 




-0.5 - 
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Fig. 17. [O/H], [N/H] and [Fe/H] vs radius for Model 7 (black 
solid line) and for Model COS (red dashed line). The data are 
the same of Figure |3] 



for the Galactic disk together with a threshold in the gas 
density for the star formation rate. In particular, the inside- 
out formation is important to reproduce the right slope 
of the abundance gradients in the inner disk, whereas the 
threshold gas density is important to reproduce the slope 
of the gradients in the outer disk. Models with a constant 
timescale for the disk formation (no inside-out) cannot re- 
produce the slope in the inner disk. On the other hand, mod- 
els with no inside-out mechanism can well reproduce the 
distributions of the star formation rate and the gas density 



along the disk. A good way of testing the inside-out disk 
formation is to test whether there is a size and luminos- 
ity evolution of disks at high redshift. Roche et al. (1998), 
by studying galaxies in the reshift range 2.5-3.0, concluded 
that the larger increase of the surface brightness relative to 
the luminosity at high redshift is an indication of an inside- 
out disk formation. They tested their results by means of 
the Chiappini et al. (1997) model. Another important is- 
sue concerns the existence of a threshold gas density for 
star formation. Recently, Boissier et al. (2006), by means 
of Galex data, suggested that the threshold might not exist 
or be much lower that thought before. Note, however, that 
the low level of star formation observed with GALEX at 
the outer edges of disks does not produce OB associations, 
hence no massive stars which are main responsible for en- 
riching the ISM. Therefore, the abundance gradients really 
measure the threshold for massive star formation. 

- In the framework of models with no inside-out mechanism 
such as the model proposed by Colavitti et al. (2008) we 
tested the effect of the efficiency of star formation vary- 
ing with the galactocentric radius, being higher in the in- 
nermost than in the outermost regions of the Galactic disk. 
This assumption, even without the threshold in the surface 
gas density, can produce gradients with the right slope both 
in the inner and outer regions of the Galactic disk but it fails 
in reproducing the gas density distribution along the disk. 

- A cosmological model without threshold, with a star for- 
mation efficiency changing with radius and with a very 
simple inside-out prescription, can well reproduce the SFR 
along the disk and the behavior of the surface gas density. 
Moreover it can well reproduce the outer abundance gradi- 
ents. In the inner part of the disk it produces too flat gradi- 
ents, even if the slope is more pronunced then that of Model 
C08. 
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- Therefore, we conclude that to reproduce at the same time 
the abundance, star formation rate and surface gas density 
gradients along the Galactic disk it is necessary to assume 
an inside-out formation for the disk. The threshold in the 
gas density is not necessary and the same effect could be 
reached by assuming a variable efficiency of star formation. 
The inclusion of radial flows without an inside-out forma- 
tion of the disk would have the same effect as the theshold 
gas density without inside-out effect (Portinari & Chiosi, 
2000). 

- It is clear from our numerical simulations that while the 
process of formation of the Galactic disk has a major im- 
pact on the formation of abundance gradients, the specific 
adopted SF law has a major effect on the gas distribution 
along the disk. This effect could be used to infer the SF law 
from the gas distributions along disks of spirals. 

- As far as the Colavitti et al. (2008) model is concerned, we 
conclude that our cosmological simulation has not a suf- 
ficient resolution to be able to capture the inside-out for- 
mation of a galaxy-sized DM halo. High resolution cos- 
mological simulations and/or resimulation of a single DM 
halo should be performed in order to understand if a suit- 
able baryon accretion law can be inferred from collisionless 
simulations, or if the full gas physics and sub-grid treat- 
ment of astrophysical processes as star formation and en- 
ergy feedback from SNe is unavoidably needed to get a re- 
aUstic description of the formation of a galaxy. 
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